Magnetic Anisotropy of Single Mn Acceptors in GaAs in an External Magnetic Field 
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We investigate the effect of an external magnetic field on the physical properties of the acceptor 
hole states associated with single Mn acceptors placed near the (110) surface of GaAs. Cross- 
sectional scanning tunneling microscopy images of the acceptor local density of states (LDOS) show 
that the strongly anisotropic hole wavefunction is not significantly affected by a magnetic field 
up to 6 T. These experimental results are supported by theoretical calculations based on a tight- 
binding model of Mn acceptors in GaAs. For Mn acceptors on the (110) surface and the subsurfaces 
immediately underneath, we find that an applied magnetic field modifies significantly the magnetic 
anisotropy landscape. However the acceptor hole wavefunction is strongly localized around the Mn 
and the LDOS is quite independent of the direction of the Mn magnetic moment. On the other 
hand, for Mn acceptors placed on deeper layers below the surface, the acceptor hole wavefunction 
is more delocalized and the corresponding LDOS is much more sensitive on the direction of the 
Mn magnetic moment. However the magnetic anisotropy energy for these magnetic impurities is 
large (up to 15 meV), and a magnetic field of 10 T can hardly change the landscape and rotate the 
direction of the Mn magnetic moment away from its easy axis. We predict that substantially larger 
magnetic fields are required to observe a significant field-dependence of the tunneling current for 
impurities located several layers below the GaAs surface. 
PACS numbers: 75.50.Pp 



14 I. INTRODUCTION 

15 Magnetic semiconductors have attracted strong atten- 

16 tion in the last decade because of their potential to 

17 combine opto-electronic and magnetic properties in spin- 
is tronic devices. The most commonly investigated material 

19 as a magnetic semiconductor is GaAs doped with tran- 

20 sition metal Mn-impurities. Mn acts as an acceptor in 

21 GaAs and its magnetic properties are mainly determined 

22 by the magnetic moment of the half filled d-shell [1]. In 

23 highly Mn doped GaAs, the observed ferromagnetism in 

24 GaMnAs has been shown to be hole mediated [2|, Q , as a 

25 result of exchange coupling between the p-like acceptor 

26 holes residing in the valence band and the electrons in the 

27 d-shell which we will refer to as the Mn core from now on. 

28 On the other hand, for applications in spintronic devices, 

29 it is important to investigate methods to read, set and 

30 manipulate the magnetic orientation of the Mn core, es- 

31 pecially at the level of a single Mn impurity. Spectacular 

32 results have been achieved with optical polarization and 

33 manipulation of low Mn doped GaAs/AlGaAsquantum 

34 wells [4] and single Mn doped quantum dots [5|, Q . Other 

35 important work in the field of single spin reading and 
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36 manipulation has been done for single nitrogen- vacancy 

37 centers in diamond [7]. 

38 In this paper, we investigate low-concentration Mn- 

39 doped GaAs. Because Mn has strongly coupled mag- 

40 netic and electric properties, spin manipulation by elec- 

41 trie fields has been suggested as a possibility in addition 

42 to manipulation by magnetic and optical fields. Cross- 

43 sectional scanning tunneling microscopy (X-STM) has 

44 been used in the past to study the Mn acceptor wave 

45 function at the atomic scale and to manipulate its charge 

46 state. The experimental study of the Mn acceptor wave- 

47 function by X-STM showed a strongly anisotropic shape 

48 of the acceptor wavefunction [8[ as was predicted by tight 

49 binding calculations [9[. These experimental and theo- 

50 retical results proved that the observed anisotropy of the 

51 acceptor wavefunction is due to the cubic symmetry of 

52 the GaAs crystal. Additional studies showed that the 

53 anisotropy of the Mn acceptor wavefunction is also in- 

54 fluenced by (local) strain due to a nearby InAs quantum 

55 dot [lO| or the relaxation of the surface [11]. 

56 These results indicate that STM can also be an excel- 

57 lent tool to investigate the effects of the orientation of the 

58 magnetic moment of the Mn core on the acceptor wave- 

59 function. In fact, theoretical work|12],|13| has predicted 

60 that the local density of states (LDOS) of the acceptor- 

61 hole wavefunction can depend strongly on the direction 



62 of the Mn moment. Since the LDOS is directly related 

63 to the tunneling current, these predictions suggest that 

64 it might be possible to control the STM electric current 

65 by manipulating the individual Mn core spin, for exam- 

66 pie with an external magnetic field. An X-STM and X- 

67 STS study of the energetic level of Mn close to the GaAs 

68 [110] cleavage surface has already shown that the 3-fold 

69 degeneracy of the J=l_CTound level is split because of 

70 the reduced symmetry [ij . Magnetic-field manipulation 

71 and control of atomic spins is presently undergoing fast 

72 progress, showing great promise to selectively address in- 

73 dividual atoms [151 . Control of atomic spin, combined 

74 with the aforementioned sensitivity of the STM current 

75 on the dopant magnetic moment direction, could be a 

76 crucial step in realizing multifunctional spin-electronic 

77 devices based on individual atoms. Apart from address- 

78 ing electrical properties of single magnetic dopants, STM 

79 has been shown to be also well capable of positioning in- 

80 dividual dopants within a semiconductor surface [16:, 17]. 

81 In this paper we will use STM to explore the effect of an 

82 external magnetic field on the magnetic orientation of the 

83 magnetic moment of a single Mn impurity in dilute Mn 

84 doped GaAs and compare the results with tight-binding 

85 model calculations. In Section [III we present a review of 

86 the theoretical work that has been published in 2 papers 

87 [l2| , [13| . These calculations are based on a tight binding 

88 model and show that a change in the spin orientation 

89 of the Mn core can sometimes give rise to a detectable 

90 change in LDOS of the Mn acceptor wavefunction. In 

91 Section imi we present experimental results of STM mea- 

92 surements on single Mn impurities in GaAs in a magnetic 

93 field. We will show that the LDOS of the Mn acceptor 

94 wavefunction is not significantly modified by magnetic 

95 fields up to 6 Tesla. In Sections |TVl and |Vl we present 

96 theoretical results of tight binding modelling of Mn in 

97 GaAs where a magnetic field has been explicitly included 

98 in the Hamiltonian. These calculations support our ex- 

99 perimental observations and show that a dependence of 
LDOS on external magnetic is in fact expected only for 
Mn acceptors placed several layers below the GaAs (110) 
and can be detected only with stronger magnetic fields 
than the ones presently available. 



104 II. REVIEW 

105 Tang et al. [l2| and Strandberg et al. [13| have re- 

106 ported results of calculations of the dependence the Mn 

107 acceptor hole wavefunction on the orientations of the Mn 

108 magnetic moment. The paper by Tang et al. [121 de- 

109 scribes the Mn LDOS in bulk GaAs with an sp^ tight 
no binding model in which the Mn core spin is taken in cal- 

111 culation by a spin dependent term in the potential at 

112 the four nearest neighbor sites in a zinc-blende crystal. 

113 It is found that the energy spectrum of the Mn is inde- 

114 pendent of the Mn core spin orientation. However, the 

115 LDOS of the Mn is found to be depending on the Mn 

116 spin orientation. A qualitative description of this depen- 



dence is given in terms of spin-orbit coupling between the 
spin of the Mn core and the orbital character of the Mn 
acceptor hole. In absence of spin-orbit interaction, the 
LDOS of the Mn acceptor state would have the same T^^ 
symmetry as the surrounding zinc-blende crystal. How- 
ever, the spin-orbit coupling is taken into account and the 
symmetry of the Mn acceptor wavefunction is reduced. 
The contour surface of the acceptor LDOS for various 
Mn core spin directions show that in general, the LDOS 
has an oblate shape with the short axis aligned with the 
Mn core spin axis. For a quantitative comparison with 
X-STM experiments, cross sectional views of the LDOS 
are calculated in the (110) plane. The largest variation 
in the cross sectional images of the LDOS is seen when 
the Mn core spin direction changes from [001] or [110] to 
[110]. A variation in LDOS of up to 90% is predicted by 
these tight binding calculations when the Mn core spin 
switches from parallel to perpendicular to the (110) sur- 
face. There is also a small difference of 15% in the LDOS 
when the Mn core spin is aligned in the two directions 
parallel to the (110) plane. When the spin of Mn core can 
be changed with an external magnetic field and possibly 
with ESR techniques [i5|,|l8|, the differences in the LDOS 
are expected to be visible in an X-STM experiment. 

This model gives a good description of Mn in bulk 
GaAs but the effect of the cleavage surface is completely 
neglected. In fact it has been shown experimentally [ll| 
that the wavefunction of a Mn near the (110) cleavage 
surface can be strongly affected by the strain from the 
surface relaxation. In the same paper, bulk tight binding 
calculations support the observation of a broken symme- 
try near the surface. The surface is taken into account by 
applying a uniform strain to the bulk model by shifting 
the Ga lattice with respect to the As lattice. The calcu- 
lation results presented in that paper are the average of 
different Mn core spin orientations. In Fig. (TJ the same 
results are presented but for individual Mn core spin ori- 
entations. Fig. [T]was unpublished in this form. A clear 
difference in LDOS can be observed when the Mn core 
spin changes its orientation from the hard axis to the 
easy axis. 

In the paper by Strandberg et al. [13] the recon- 
structed surface is taken into account for the calculation 
of the LDOS dependence of the Mn acceptor state on 
the Mn core spin orientation using a tight binding model 
which includes the exchange interaction, spin-orbit cou- 
pling and Coulomb interaction. This makes a direct com- 
parison with X-STM experiments more justified and the 
results indeed show the same experimentally observed 
breaking of the symmetry of the wavefunction due to the 
near presence of the surface. In Ref. 13 Mn acceptors in 
bulk GaAs (neglecting the surface) have also been con- 
sidered. For Mn in bulk GaAs, the energy level of the 
Mn state calculated for different orientations of the Mn 
core spin shows a small magnetic anisotropy, in contrast 
to the results of Ref. [l2|, where no magnetic anisotropy 
was found for Mn in bulk. The easy axis in Ref. T3 for 
the Mn core spin is oriented along the [001] direction 
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FIG. 1. (a) Magnetic Anisotropy Energy (MAE) of a Mn 
acceptor according to a tight binding calculation for strained 
bulk material. The energy level difference between the easy 
and the hard axis is about 23 meV, based on a uniform strain 
estimated in Ref. [ir]. The angles 6 and (j) have the same 
definition as in Fig. |4l (b) Mn LDOS at five atomic layers 
from the Mn position when the core spin is oriented along the 
(b) easy axis or (c) hard axis. 



whereas the hard axis is found to be lying in the (001) 
plane. The energy barrier between the hard axis and the 
easy axis is found to be 4.35 meV which is very small 
in comparison with the Mn binding energy in GaAs (113 
meV). At first, the presence of a magnetic anisotropy is 
surprising since there is no difference between the [001] 
and [010] or [100] directions in a Zinc-Blende crystal. 
The observed anisotropy can be explained by the use of 
periodic boundary conditions on finite clusters used in 
this paper [13| . The influence of other Mn atoms in the 
area may indeed introduce a small magnetic anisotropy 
and thus the observed magnetic anisotropy of Mn in bulk 
GaAs is artificial. Indeed, more recent calculations car- 
ried on out on much larger clusters show that the bulk 
magnetic anisotropy decreases monotonically with clus- 
ter size, down to a fraction of a meV for the largest clus- 
ters of 40,000 atoms 0- 

On the other hand, the calculation of the LDOS for 
Mn in bulk GaAs in Ref. [13il shows good similarity with 
the calculations in Ref. [12|. The LDOS is found to be 
spreading in the direction perpendicular to the Mn core 
spin axis. The change in the shape of the LDOS is ex- 
plained in terms of the p^^, p^y and P;^ character of the 
Mn acceptor hole. For different orientations of the Mn 
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core spin, different components in the character domi- 
nate. When the Mn core spin direction is changed from 
[110] to [110] a drop in LDOS of 74% is observed at 4 
atomic layers from the Mn position. This drop in LDOS 
is 25% when the core spin direction changes from [110] 
to [001], which is again in good agreement with the other 
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calculations in 

In Ref. [l3| similar calculations have been done for Mn 
in or below the GaAs (110) surface layer. For Mn at the 
surface and the first subsurface layer, a strong localiza- 
tion of the LDOS is observed and a magnetic easy axis in 
the [111] direction is found. The difference in LDOS for 
different Mn core spin orientations is negligible. Thus in 
an X-STM experiment, we expect to see no effect of the 
magnetic field on the Mn atoms very close to the surface. 

For Mn atoms deeper below the (110) surface, 
the LDOS becomes more extended and the magnetic 
anisotropy shows a complex behavior for subsequent 
depths. However, from the fourth layer beneath the (110) 
surface and deeper, one can recognize the emergence of 
an easy plane with its normal in the [110] direction. The 
anisotropy energy is found to be at least 15 meV. Images 
of the (110) surface LDOS show that there is an increas- 
ing difference in LDOS for an increasing depth when the 
Mn core spin changes from the easy axis to the hard 
axis. For Mn atoms placed on fourth subsurface layers 
and deeper, the difference in LDOS varies between 40% 
and 82%. 

In summary, both Refs. Il2| and |l3| have treated the 
behavior of the Mn acceptor hole LDOS in the (110) 
plane for different Mn core spin orientations. In both 
papers it is found that when the Mn core spin direction 
is changed from [110] to [110], a drastic change in the 
LDOS is taking place. The inclusion of the cleavage sur- 
face relaxation has resulted in similar observations. 

The mechanism for the magnetic anisotropy in Refs.[l2l 
and Uj is the same — the presence of the surface, or 
strain, lowers the energy of an orbital wave function with 
quantization axis along a specific direction, and the spin- 
orbit interaction (which correlates the spin axis with the 
orbital axis) causes that preferred orbital direction to 
select a preferred spin axis. The effective energy asso- 
ciated with the correlation between spin axis and or- 
bital axis is of the same order as the binding energy 
(Refs. H, m, and [H found it to be - 40 meV). At mag- 
netic fields required to overcome the magnetic anisotropy 
energy the magnetic length is of order 3 nm, which is 
three times larger than the effective Bohr radius of the 
acceptor [r^ 1 nm). Therefore the overall distortion of 
the acceptor state wave function due to the direct effect 
of the magnetic field on the orbital wave function is small 
compared with the spin-orbit term. What is not certain, 
however, is whether the effect of the magnetic field on the 
acceptor state wave function can substantially change the 
magnetic anisotropy; this will be examined in Section IVl 

In an X-STM experiment, one can also check the re- 
sults of these calculations by applying a magnetic field 
perpendicular and parallel to the (110) cleavage plane 



258 and by measuring the Mn contrast, which can change 

259 with a factor as high as 90%. In the next section, we dis- 

260 cuss the X-STM experiments that have been performed 

261 to observe the predicted effects. 



262 III. EXPERIMENTS 

263 A Mn doped layer of 500 nm thick was grown on a 

264 p-type GaAs substrate at high temperature with Mn 

265 concentrations of about 3x10^^ cm~^. The experiments 

266 are performed with an Omicron Cryogenic STM oper- 

267 ating at a base temperature of 2.5 K. A magnetic field 

268 vector can be applied with fields of up to 6 T in the 

269 z-direction only or max. 2 T in the z-direction together 

270 with max. 1 T in the x- and y-directions. The x-, y- and 

271 z-direction are indicated in figure [2^ where 2 Mn atoms 

272 at different depths below the cleaved surface are visible. 
27a The magnetic field is indicated in the vector notation in 

275 units of T: B={B^,By,B^). 

276 From Ref. [11| , we estimate that Mn A is approximately 

277 8 atomic layers below the cleavage surface and that Mn 

278 B is at about 5 atomic layers below the cleavage surface. 

279 In [l3|, a change in contrast of 40% is predicted for a 

280 Mn A at 8 layers below the cleavage surface when the 

281 Mn core spin changes from the [110] direction to the 

282 [110] direction. For Mn B at 5 atomic layers beneath 

283 the cleavage surface, a change of 60% is predicted 

284 when the Mn core spin direction changes from the [110] 

285 direction to the [110] direction. As can be seen from 

286 the comparison of figures [2^ and [2]3, there is no change 

287 at all in the Mn contrast for both Mn atoms when the 

288 magnetic field is changed from 1 T in the y direction to 

289 -6 T in the z-direction. In figure [3l a more quantitative 

291 comparison is made by looking at the contrast of the Mn 

292 atoms in different magnetic fields through the dashed 

293 lines in figure [SJ). Also in these plots, it can be seen that 

294 for both Mn atoms, there is no difference at all in the 

295 contrast for different orientations of the magnetic field. 

296 For Mn B, the plots for B=(0,1,0) and B=(0,0,-6) are 

297 slightly different from the rest because of a small tip 

298 modification that has taken place. The tip has become 

299 slightly less sharp in the scan direction (the [001] 

300 x-direction) and this difference is noticed when sharper 

301 objects like Mn B are imaged. Mn A has FWHM of 

302 about 4.5 nm in the scan direction, while Mn B has a 

303 sharper feature with a FWHM of about 2.0 nm. The 

304 different FWHM of the Mn features has been related to 

305 the depth below the GaAs surface [ll| 



IV. THEORETICAL MODEL 

We model theoretically substitutional Mn impuri- 
ties in GaAs following the procedure put forward in 
Ref. [Isl. Our second-quantized tight-binding Hamilto- 




FIG. 2. 13x13 nm X-STM images of two Mn atoms at differ- 
ent depths below the (110) cleavage plane. The x direction 
corresponds with the crystallographic [001] direction and the 
y-direction corresponds with the [110] direction. The images 
are taken at +1.4 V and 50 pA at a temperature of about 
2.5 K. In a) a magnetic field of 1 T is oriented in the [110] 
y-direction and in b) a magnetic field of -6 T is oriented in 
the [110] z-direction. 



nian for (Ga,Mn)As takes the following form: 
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308 where i and j are atomic indices that run over all atoms, 

309 m runs over the Mn, and n[m] over the nearest neighbors 

310 of Mn atom m. ji and v are orbital indices and a is a 

311 spin index. The first term in Eq. ([1]) contains the near- 

312 neighbor Slater-Koster tight-binding parameter s 1201 121| 

313 that reproduce the band structure of bulk GaAsj22| and 
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FIG. 3. a) Contrast of Mn A along the [001] direction indi- 
cated with a dashed hne in figure [2]3. b) same plot for Mn B. 
Mn A has FWHM of about 4.5 nm and Mn B has a FWHM 
of about 2.0 nm 



314 that are rescaledp^ |23|, |2J when needed to account for 

315 the buckling of the (110) surface. 

316 The second term implements the antiferromagnetic ex- 

317 change coupling between the Mn spin Clrn (treated as 

318 a classical vector) and the nearest neighbor As p-spins 
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The exchange coupling 



319 Sn 

320 Jpd = 1.5 eV has been inferred from theory [25| and ex- 

321 periment ^] . As a result of this term the acceptor hole 

322 that is weakly bound to the Mn will become spin polar- 

323 ized. This model contains only s and p orbitals, and the 

324 effect of the Mn 3d^ electrons is encoded in the exchange 

325 term. 

326 Next, we include an on-site spin-orbit one-body term, 

327 where the renormalized spin-orbit splittings are taken 

328 from Ref. [22|. Spin-orbit coupling will cause the total 

329 energy to depend on the Mn spin direction, defined by a 

330 collinear variation of Clm • 

331 The fourth term is a long-range repulsive Coulomb part 



332 that is dielectrically screened by the host material. To 

333 account in a simple way for weaker dielectric screening 

334 at the surface, the dielectric constant e^ for a Mn on 

335 the surface is reduced from the bulk GaAs value 12 to 

336 6 for the affected surface atoms. This crude choice is 

337 qualitatively supported by experimental results [13, 123 • 

338 The last term Vcorr is a one-particle correction poten- 

339 tial for the Mn central cell. This term is the least known 

340 and understood theoretically. It consists of on- and off- 

341 site parts, Korr = V^n + V^ff which influence the Mn ion 

342 and its As nearest neighbors respectively. The on-site 

343 Coulomb correction is estimated to be 1.0 eV from the 

344 ionization energy of Mn. The off-site Coulomb correc- 

345 tion affects all the nearest-neighbor As atoms surround- 

346 ing the Mn ion and together with the exchange interac- 

347 tion, it reflects primarily the p-d hybridization physics 

348 and is the parameter that in the model primarily con- 

349 trols the binding energy of the hole acceptor state. The 

350 off-site Coulomb correction value is set by tuning the po- 

351 sition of the Mn-induced acceptor level in the bulk to 

352 the experimentally observed position |28l-l31| at 113 meV 

353 above the first valence band level. The value thus ob- 

354 tained is Kff = 2.4 eV. When the Mn impurity is on the 

355 GaAs surface, the value of V^ff is reduced to ensure that 

356 the position of the acceptor level is consistent with the 

357 value attained via STM spectroscopy. 

358 The off-site Coulomb correction is in fact a repulsive 

359 potential for the electrons. If we use the bulk value (2.4 

360 eV) for the surface, the acceptor level lies deep in the gap 

361 at 1.3 eV above the valence band, which means the accep- 

362 tor wave function is now much more localized around the 

363 Mn than its bulk counterpart. In order to guarantee the 

364 experimentally observed position for the acceptor level, 

365 0.85 eV [161, we have to decrease this repulsive potential 

366 for the electrons, which causes the hole wave function 

367 to be less localized with a corresponding smaller binding 

368 energy. 

369 The electronic structure of GaAs with a single substi- 

370 tutional Mn atom is obtained by performing a super-cell 

371 type calculation with a cubic cluster of a few thousands 

372 atoms and periodic boundary conditions in either 2 or 3 

373 dimensions, depending on whether we are studying the 

374 (110) surface or a bulk-like system. The (110) surface 

375 of GaAs is simplified from both theoretical and experi- 

376 mental points of view, by the absence of large surface re- 

377 construction. In order to remove artificial dangling-bond 

378 states that would otherwise appear in the band gap, we 

379 include relaxation of surface layer positions following a 

380 procedure introduced in Refs. j23l and |2^ . For more de- 

381 tails the reader is referred to Ref. |13| . 

382 We would like to emphasize that the strength of the 

383 off-site Coulomb correction is the only important fitting 

384 parameter of the model, and its value is fixed once for 

385 all by the procedure described above. All the other pa- 

386 rameters in Eq. [T] are either determined by theoretical 

387 considerations, or for the cases when this is not possible 

388 (e.g. short- range onsite potential) their values are ex- 

389 tr acted from experiment. In any case, they affect weakly 
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the properties of the acceptor level. Once the parame- 
ters of the Hamiltonian of Eq. [1] are chosen in the way 
indicated above, the model has to be viewed as a micro- 
scopic description, with predictive power, of the proper- 
ties of Mn impurities in GaAs surfaces and subsurfaces. 
In this sense the model of Ref. [isl has been quite suc- 
cessful in capturing some of the salient features of the 
STM experiments ^4, |32j, probing the Mn-dopant ac- 
ceptor hole near the GaAs (110) surface. For example, it 
correctly describes the dependence of the acceptor bind- 
ing energy [1^ and the shape of the hole wave function 
|32| on the layer depth below the surface on which the 
magnetic dopant is positioned. The model also makes a 
prediction on how the magnetic anisotropy barrier for the 
Mn- impurity-hole magnetic complex changes as a func- 
tion of the layer depth. These predictions can be indi- 
rectly checked by the magnetic-field studies that are the 
main scope of the present paper. 

In order to study the response of the system to an 
external magnetic field, we introduce the Zeeman term 
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410 where the first term runs over all s and p orbitals of 

411 all atoms, and the second term represents the coupling 

412 of the magnetic field with the magnetic moment of the 

413 Mn impurities, treated as a classical vector. Here jib = 

414 ^ = 5.788 X 10"^ meV T"^ is the Bohr magneton, Qs = 

415 2, and we follow the incorrect but common convention 

416 that spins and magnetic moments are parallel to each 

417 other [33| . Therefore in the paper we will loosely refer to 

418 the direction of Cl as the direction of the Mn magnetic 

419 moment. 



V. THEORETICAL RESULTS AND 
DISCUSSION 
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We start by analyzing the magnetic anisotropy prop- 
erties for one Mn at the (110) GaAs surface layer and 
the immediate subsurface layers, and see how these are 
modified by the presence of an external magnetic field 
of a few Tesla. The magnetic anisotropy landscape as a 
function of Cl for one Mn at the surface and the first 9 
subsurfaces has been studied in detail in Ref. [l3|. Typ- 
ically the system has an uniaxial anisotropy with two 
minima separated by an energy barrier. We will refer to 
the Cl direction of minimum energy as the easy direction 
and the one of maximum energy as the hard direction. 

We first consider the case of one Mn impurity at the 
(110) surface. To facilitate the comparison with the 
case in which a magnetic field is present, we recalcu- 
lated and plotted here anisotropy landscapes and LDOS 
in the absence of the magnetic field, originally published 
in Ref. |l3|, using an improved code. 
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FIG. 5. Color online - The magnetic anisotropy energy 
(MAE) landscape and the Mn- acceptor- level local density of 
states (LDOS) for one Mn at the [110] surface, (a) The MAE 
in the absence of an external magnetic field, as a function 
of the angles 6 and (f) defining the direction of the Mn spin 
Cl. The barrier (hard) direction is marked with a circle and 
the minimum energy (easy direction) with a square, (b) and 
(c) LDOS of the Mn acceptor level when the Mn magnetic 
moment points in the easy and hard direction respectively, 
as defined in panel (a), (d) The MAE in the presence of a 
6 T external magnetic field applied along the hard direction 
{0 = 37r/4, = 7r/4). (e) The LDOS in the presence of a 6 T 
magnetic field. Here the Mn magnetic moment is along the 
easy direction determined by the landscape (d) modified by 
the presence of the external field. The barrier (hard) direc- 
tion is marked with a circle and the minimum energy (easy 
direction) with a solid line. 



439 In Fig. [5ja) we plot the anisotropy energy landscape 

440 in the absence of the magnetic field, as a function of the 

441 angles and (f) defining the direction of Cl. The coordinate 

442 system used for this and the other plots in the paper has 

443 = parallel to the [001] axis, {0 = 7r/2, = 0) parallel 

444 to [100], and (l9 = 7r/2, = 7r/2) parallel to [010]. See 
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Fig.H 

The anisotropy landscape displays two minima, identi- 
fying the easy direction [111], separated by an energy bar- 
rier of the order of 1 meV. Note that these tight-binding 
results of the magnetic anisotropy of a Mn at the (110) 
GaAs surface are consistent with recent first-principles 
estimates [34|. Panels (b) and (c) of Fig. [5] show the LDOS 
for the Mn acceptor state when the Mn spins point along 
the easy and hard direction respectively, determined from 
the landscape in (a). As discussed in Sec. II and shown 
clearly in the figures, the acceptor state wavefunction for 
a Mn on the surface is very localized around the impurity, 
and the dependence of the LDOS on the Mn spin orien- 
tation is negligible. The acceptor wavefunction, itself 
strongly anisotropic, seems to be completely decoupled 
from the orientation of the Mn magnetic moment. 

Fig. [5fd) and (e) show the effect of a 6 T magnetic 
field on the anisotropy and LDOS respectively, when the 
field is applied in the hard direction of the anisotropy 
landscape in (a). 

We can see that the magnetic field changes consider- 
ably the anisotropy landscape, which has now an easy 
axis at ^ = TT (the direction of the field). Note that in 
the presence of the field the anisotropy barrier has in- 
creased up to 4 meV. The LDOS in Fig. [5fe) is now 
calculated for (l pointing along the new easy axis, deter- 
mined by the magnetic field. Despite the strong change 
in the anisotropy landscape brought about by the mag- 
netic field, the acceptor LDOS is essentially identical to 
the one calculated in the absence of the field, in agree- 
ment with the experimental results. 

Before continuing our LDOS analysis, it is useful to 
consider how the anisotropy-energy barrier depends on 
the Mn-impurity depth from the (110) surface. In Fig. [6] 
we plot the largest value of the anisotropy energy barrier 
as a function of the subsurface layer index (layer is the 
(110) surface). In general, in the absence of a magnetic 
field (red dots in the picture) the anisotropy barrier in- 
creases with Mn depth, reaching a maximum of 15 meV 
for layers 4 and 5. It then starts to decrease and it should 
eventually reach a very small value corresponding to the 
case where the Mn is effectively in the bulk. For the finite 
clusters that we have considering here (20 layers in the 
z-direction), the anisotropy remains large also when the 
impurity is effectively in the middle of the cluster (corre- 
sponding to layer 9 from the surface). Bulk calculations 
on considerably larger clusters show that the anisotropy 
for impurities in the middle of the clusters does decrease 
to a fraction of one meV|19l|. For these larger clusters 
the magnetic anisotropy of the Mn positioned in on lay- 
ers > 8 is expected to decrease a bit with cluster size. 
However the qualitative behavior of the first 7-8 layers 
shown in Fig. [6l and the corresponding numerical values 
of the magnetic anisotropy are controlled by the vicinity 
to the surface and as such should not depend strongly on 
cluster sizeflQj. 

Layer 1 (the first subsurface layer) is a special case in 
the sense that the anisotropy is very small, on the or- 
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FIG. 6. Color online - The maximum MAE barrier height as 
a function of the Mn depth. Red dots are the MAE barrier 
height in the absence of an external magnetic field, while blue 
dots represent the height in the presence of a 6 T external 
magnetic field. 



der of 0.1 meV. The first subsurface represents the cross 
over from the case in which the Mn is at the surface, 
with three nearest neighbor As, to a bulk-like environ- 
ment characterized by four nearest neighbor As atoms. 
The properties of the acceptor level found in STM exper- 
iments for a Mn positioned on this subsurface are also 
quite anomalous |35|. When a magnetic field of 6 T is 
applied along the hard direction (blue dots in Fig. [6]) 
the anisotropy barrier increases by a couple of meV. The 
exception is again the first subsurface (layer 1), whose 
anisotropy is now completely controlled by the magnetic 
field and behaves in a similar way to the surface layer. 
The behavior of the first subsurface anisotropy landscape 
is shown explicitly in Fig. [7] (a), (d). 

As for the case of a Mn atom placed at the (110) sur- 
face, the acceptor LDOS for a Mn on the first subsurface 
(see Fig.[71^b), (c) ) is completely insensitive to the direc- 
tion of the Mn magnetic moment. Again a 6 T magnetic 
field, which is able to completely modify the magnetic 
anisotropy landscape and orient the Mn moment parallel 
to its direction, does not have any detectable effect on 
the acceptor wave function, as shown in Fig. [Tl^d). 

As the Mn is placed in successively deeper layers below 
the surface and the acceptor wavefunction becomes less 
localized around the impurity, the situation changes. In 
Figs. [8] and [9] we plot the anisotropy landscape and the 
acceptor LDOS for the fourth and fifth subsurface (layer 
4 and 5 below the surface) respectively. As we discussed 
before, when the direction of the Mn moment is forced 
to point in the hard direction (panel (c) of Fig. [8j) the 
LDOS around the Mn increases sensibly. The two cases, 
easy (panel b) and hard axis LDOS are now clearly dis- 
tinguishable. Since the acceptor wavefunction is always 
normalized, an increase of the LDOS in the core region 
implies that the acceptor wavefunction is considerably 
more localized when the Mn magnetic moment points in 
the hard direction. On the other hand, in contrast to 
the surface and the first subsurface, the energy barrier in 
this case is considerably larger. A magnetic field of the 
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FIG. 7. The MAE landscape and the Mn- acceptor- level 
LDOS for one Mn in the first subsurface (i.e., one layer below 
the [110] surface), (a) The MAE in the absence of an exter- 
nal magnetic field. The barriers (hard) directions are marked 
with a circle, (b) and (c) the Mn acceptor LDOS for the case 
in which the Mn magnetic moment points in the easy and 
hard direction respectively, (d) The MAE in the presence of 
a 6 T external magnetic field pointing along the (original) 
hard direction {0 — O^cj) — n). The minimum energy (easy 
direction) is shown with a solid line, (e) The Mn acceptor 
LDOS in the presence of a 6 T magnetic field. Here the Mn 
magnetic moment points in the new easy direction determined 
by the magnetic field, as shown in (d). The colorscale in (b) 
and (c) is the same as in (e). 




FIG. 8. The magnetic anisotropy energy (MAE) and the 
LDOS for one Mn in the fourth subsurface (fourth layer below 
the [110] surface), (a) The MAE in the absence of an external 
magnetic field, (b) and (c) the LDOS of Mn acceptor level 
for the case that Mn magnetic moment points in the easy and 
hard direction respectively, (d) The MAE in the presence of 
a 6 T external magnetic field which points along the hard 
direction(6' = 7r/2, = 37r/4). (e) The LDOS in the case that 
magnetic moment points in the easy direction in the presence 
of a 6 T magnetic field. The barrier direction in (a) and (d) 
is marked with a circle and the easy direction with a square. 
The colorscale in (b) and (c) is the same as in (e). 



order of those applied experimentally are now not strong 
enough to modify appreciably the anisotropy landscape. 
This can be seen by comparing panel (a) - no magnetic 
field - with panel (d), where magnetic field of 6 T is ap- 
plied in the hard direction. Consequently, the direction of 
the easy axis is only slightly modified in the presence of a 
magnetic field, and as a result, the corresponding accep- 
tor LDOS appears now very similar to the zero-magnetic 
field case [panel (b)]. This is again in agreement with the 
experiments presented in this paper. 

As mentioned before, the experiments presented in [l^ 
showed the energy level splitting of Mn in GaAs close to 
the cleavage surface. For a typical Mn position at 5th 
subsurface layer, a total splitting of 14 meV is found be- 
tween the 3 peaks which are attributed to the different 
projections of the total momentum J=l which is the re- 
sult of anti-ferromagnetic coupling between the 5/2 Mn 
core spin and 3/2 Mn acceptor total angular momentum. 
In Fig. [9l it can be seen that the MAE is indeed about 
15 meV which corresponds well with the findings in [l^ . 

In Fig. [U the MAE as calculated in another tight bind- 
ing calculation (strained bulk GaAs) is about 23 meV. 
This is more than the 15 meV of the supercell calcu- 
lations (Fig. [9|) possibly because of the overestimated 



566 strain or its assumed uniformity. 

567 We conclude that, although the LDOS of deep- 

568 subsurface Mn acceptors is in principle strongly depen- 

569 dent on the Mn magnetic moment direction, its actual 

570 manipulation with an external magnetic field is not suit- 

571 able at field strengths presently used in experiment. 



VI. CONCLUSIONS 



In conclusion, this work is the first systematic study 
of the effect of an external applied magnetic field on the 
acceptor properties of individual Mn impurities in GaAs. 
Specifically, we have investigated theoretically and ex- 
perimentally the effect of an external magnetic field on 
the acceptor hole wavefunction and LDOS of Mn impuri- 
ties placed near the (110) surface of GaAs. The acceptor 
LDOS is directly accessible via X-STM experiments. 

The motivation of this study was in part provided 
by previous theoretical studies which predicted that the 
LDOS in some cases strongly depends on the orientation 
of the magnetic impurity magnetic moment. The theoret- 
ically model used in this analysis is essentially parameter- 
free, once the energy of the surface acceptor state is fixed 
to reproduce the experimental value. 
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FIG. 9. As in Fig. [51 but for the fifth subsurface (fifth layer 
below the [110] surface). The colorscale in (b) and (c) is the 
same as in (e). 



588 Experimentally we find that there is no detectable dif- 

589 ference in the STM images of the acceptor hole LDOS 

590 when a magnetic field up to 6 T is applied in several 

591 directions with respect to the crystal structure. To rec- 

592 oncile theory and experiment we have carried out a the- 

593 oretical analysis of the magnetic anisotropy energy and 

594 acceptor hole wavefunction in the presence of a magnetic 

595 field. We have shown that for Mn impurities placed in 

596 deep sub- layers below the surface, the calculated mag- 

597 netic anisotropy landscape is characterized by energy 

598 barriers of the order of 10-20 meV, which are only mini- 

599 mally affected by magnetic fields used in experiment. We 

600 estimate that one needs to employ much stronger fields 

601 (on the order of tens of Tesla) to modify significantly the 

602 anisotropy landscape and rotate the magnetic moment 

603 of the impurity. This estimate is based on the idea of 

604 manipulating a spin=5/2 object with g-factor=2 with an 

605 external field to overcome an energy barrier of 15 meV. 

606 For impurities placed near the surface, the magnetic 



607 anisotropy is small enough to be considerably affected by 

608 a magnetic field of a few T. However, for this case the ac- 

609 ceptor hole LDOS is much less sensitive to the orientation 

610 of the Mn magnetic moment. The combination of these 

611 two facts seem to explain the experimental finding that 

612 the the STM images of the acceptor hole wavefunction is 

613 essentially unaffected by an external magnetic field. 

614 Our studies show that the Mn-dopant behavior close 

615 to the GaAs surface depends on the layer depth in a 

616 complex and highly non trivial way. These studies also 

617 suggest that it could be interesting to carry out a sim- 

618 ilar investigation for other magnetic dopants and other 

619 semiconductors. It might be possible that for some of 

620 these systems the acceptor wavefunction for a dopant 

621 near the surface be more delocalized and amenable to an 

622 easier manipulation by a static magnetic field, display- 

623 ing the effects originally predicted for Mn in GaAs. It 

624 should also be possible to use resonant techniques, such 

625 as those commonly used in electron spin resonance and 

626 ferromagnetic resonance, to map out the anisotropy land- 

627 scape presented here for Mn near the GaAs surface. Fi- 

628 nally, excitations of the spin that would correspond to the 

629 quantized spin in the anisotropy landscape here should 

630 be visible in inelastic tunneling spectroscopy. Thus these 

631 new predictions do not mean that Mn spin dynamics is 

632 impossible to see near the surface of GaAs, merely that 

633 it is more challenging to observe. 
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